function Gaussian_elimination(A,b)
X = horzcat(A,b);
Y = abs(X);
for i = 1:length(A)-1
    for j = i+1:length(A)
        X(j,:) = X(j,:) - X(j,i)/X(i,i)*X(i,:);
    end
end
for j = length(A)-1:-1:1
    for i = length(A):-1:j+1
        X(j,:) = X(j,:) - X(j,i)/X(i,i)*X(i,:);
    end
end
for i = 1:length(A)
    X(i,:) = X(i,:)/X(i,i);
end
X = X(:,length(A)+1);
disp(X);
end